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ABSTRACT 

Context. On 4 July 2005 at 05:52 UT, the impactor of NASA's Deep Impact (DI) mission crashed into comet 9P/Tempel 1 with a 
velocity of about lOkms"'. The material ejected by the impact expanded into the normal coma, produced by ordinary cometary 
activity. 

Aims. Based on visible and near-IR observations, the characteristics and the evolution with time of the cloud of solid particles released 
by the impact is studied in order to gain insight into the composition of the nucleus of the comet. An analysis of solid particles in the 
coma not related to the impact was also performed. 

Methods. The characteristics of the non-impact coma and cloud produced by the impact were studied by observations in the visible 
wavelengths and in the near-IR. The scattering characteristics of the "normal" coma of solid particles were studied by comparing 
images in various spectral regions, from the UV to the near-IR. For each filter, an image of the "normal" coma was then subtracted 
from images obtained in the period after the impact, revealing the contribution of the particles released by the impact. 
Results. For the non-impact coma the Afp, a proxy of the dust production, has been measured in various spectral regions. The presence 
of sublimating grains has been detected. Their lifetime was found to be ~ 11 hours. Regarding the cloud produced by the impact, the 
total geometric cross section multiplied by the albedo, SA, was measured as a function of the color and time. The projected velocity 
appeared to obey a Gaussian distribution with the average velocity of the order of 115ms"'. By comparing the observations taken 
about 3 hours after the impact, we have found a strong decrease in the cross section in J filter, while that in Ks remained almost 
constant. This is interpreted as the result of sublimation of grains dominated by particles of sizes of the order of some microns. 

Key words. Comet 9P/Tempel I - Deep Impact event - organic grains - dust ejecta cloud 



1. Introduction 

The Deep Impact mission (hereafter DI) to the Jupiter family 
comet 9P/Tempel 1 (hereafter 9P) was aimed at studying the 
cratering physics in minor bodies in the solar system and the 
primordial material preserved inside cometary nuclei. On July 4, 
2005 the impactor of the DI experiment produced a high-speed 
(about 10km s"') impact in the nucleus of 9P excavating a con- 
siderable amount of cometary material that was observed and 
measured both in-situ by the DI fly-by spacecraft and remotely 
by Earth-based instrumentation. First results of the mission are 
described in |A'Heam et al. | ( j2005l ) and [Sunshine et al. | ( |2006| l. 
Early earth-based and other space-based measurements of the 
event have been published by |Meech et al. | ( [2005[ ), [Sugita et al. | 
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(|2005] l ,[Hger et al. | ( |2005| ), |Lisse et al.1 ( |2006l l, and |Schleicherl 
|etar| ( |2005l i; 

At the European Southern Observatory (ESO) DI received 
considerable observing time allocated to observe the event at 
their Chilean observatory, at Cerro La Silla and at Cerro Paranal 
sites ( |Kaeufl et al. ||2005 a1). Here, we summarize results from the 
visible and near-IR measurements of the dust in the cometary 
coma obtained both shortly before and after the DI event. We 
focus on the dust ejecta properties such as scattering proper- 
ties, projected velocity, and spatial distribution and their evo- 
lution with time. Complementary data from the ESO DI cam- 
paign on polarimetric and mid-IR observations as well as on 
the cometary gas emission and the large-scale coma activity of 
the comet are described elsewhere (see e. g., Boehnhardt et al., 
2007). Pre-impact monitoring of the cometary activity is de- 
scribed by |Kaeufl et ar| ( |2005b| ) and [Lara et al1 ( |2006| ). 
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2. Observations 

2.1. Telescopes, instruments, filters 

The majority of the observations, described here, were per- 
formed at the European Southern Observatory (ESO) in La 
Silla/Chile using the 3.5 m New Technology Telescope (NTT) 
by switching between two focal plane instruments: EMMI (ESO 
Multi-Mode Instrument), for the visible spectral region, and 
SOFI (Son of ISAAC), for the near-IR {JHK). Both instruments 
are of focal reducer type for imaging and spectroscopic observa- 
tions. EMMI provides a field of view of 9.1 x 9.9 arcmin with a 
two-detector array in the red (400 - 1 000 nm) and of 6.2 x 6.2 
arcmin with a single detector in the blue arm (300-500nm) at 
0.32 and 0.37 arcsec/pixel resolution (using the 2x2 and 1x1 
binning options), respectively. In its large field option used for 
these observations, SOFI has a single detector of 4.9 x 4.9 ar- 
cmin field of view at 0.288 arcsec/pixel resolution. In the visible, 
narrow band filters, with bandpasses within selected wavelength 
regions of interest for cometary science, were used. In particu- 
lar, for the study of the cometary dust, the following filters with 
no or negligible gas emission in their passband were used: one 
in the ultraviolet (Uc), one in the blue (B^) and one in the red 
(Rc) spectral region. The near-IR observations were performed 
with the regular J, H, and broad band filters since in this re- 
gion the gas contamination is negligible. Table [T] gives the log 
of observations together with the list of filters used including 
the respective central wavelength and full width at half maxi- 
mum (FWHM) of the wavelength passband. Technical informa- 
tion on the La Silla telescope and instruments can be found at 
pittp ://w w w.ls .e so .org/lasilla/sciops I 

Since La Silla was clouded over during night 5-6 July 
2006, near-IR imaging of the comet was shifted to the on-going 
DI campaign at the Cerro Paranal Observatory using ISAAC 
(Infrared Spectrometer And Array Camera) at the 8.2 m unit tele- 
scope Antu of the ESO's Very Large Telescope (VLT). Due to 
the shortage of time at the end of the nightly visibility window, 
only part of the J, H and filter imaging sequence was per- 
formed. ISAAC is a focal reducer-type instrument providing a 
field of view of 2.5 x 2.5 arcmin at a pixel resolution of 0.148 
arcsec. Technical information on the VLT and the ISAAC in- 
strument can be found at http://www.eso.org/parana l/sciopsi 

2.2. Calibrations on the sky 

For calibration purposes, some photometric standards were also 
observed before and after the comet observations on each clear 
night. In the near-IR, the normal JHK^ photometric standards 
were used ( [Persson et al. [ |1998| l, while for the calibration of 
the narrow band filters in the visible, well known spectropho- 

11994) were 



tometric standards from the list by |Hamuy et al. 



measured. The required calibration frames (bias and sky flatfield 
exposures for the visible imaging and screen and/or lamp flat- 
fields with lamp illumination on and off for the near-IR) were 
obtained during daytime and/or twilight periods. 

2.3. Observing techniques 

The comet imaging was performed with the telescope tracking at 
the speed of the moving target. Jitter offsets of small amplitude 
(order of 10-30 arcsec) were applied between individual expo- 
sures through a single filter As usual for extended objects, the 
observations of the comet in the near-IR spectral region were 
interlaced by observations of the sky at an offset of ^ 8' in a 



different region of the sky. A sequence of 5 comet and 5 sky 
images were usually taken in each near-IR filter. The jitter se- 
quence typically lasted for 11-12 minutes per filter. Due to the 
mentioned shortage of time, observations on night July 5-6 with 
ISAAC consisted of only 2 comet and 2 sky images per filter. 
Observations in the visible were also repeated 5 times for each 
filter, offsetting the telescope by 10-30 arcsec. Calibration obser- 
vations (standards, sky flatfields) were performed with the tele- 
scope tracking at the sidereal rate. Daytime calibration images 
(bias, dome flats) used fixed telescope pointing. 

Since EMMI and SOFI focal plane instruments were 
mounted on the two Nasmyth foci of the NTT telescope, fully 
simultaneous observations in the near-IR and visible were not 
possible. However, the switching time between the two instru- 
ments was short (less than 15 min) and allowed us to use both 
instruments sequentially during the nightly visibility window of 
the comet. 

The summary log of observations is given in Table[T] During 
the observing period the Sun (r/,) and Earth (A) distances of the 
comet were 1.51 AU and 0.89-0.91 AU, respectively. The phase 
(Sun-Comet-Observer) angle was ^ 41°, and the position angle 
of the Sun projected on the sky at the position of the comet, was 
^ 290°. 



3. Data reduction 

3. 1 . Frame pre-processing 

For the visible imaging, all comet and standard star images 
were corrected for the bias and the flatfield. Both bias and 
flatfield maps were computed as the average of a series of 
bias and sky flatfield exposures taken during the observing in- 
terval and through the corresponding filters (for the flatfield). 
Subsequently, first-order sky background correction was applied 
by subtracting the average sky flux value measured at the four 
edges of the individual flatfielded images. 

For the near-IR data the flatfield maps were computed from 
the screen flat images in each filter with the lamp illumination on 
and off. Then, for each sequence, a median average sky+bias was 
computed from the sequence of five sky observations. The comet 
images were then reduced by subtracting the median averaged 
sky frame and by dividing the result by the flatfield for the cor- 
responding filter. Finally, comet images for each filter/sequence 
in the visible and near-IR were obtained as the median average 
of the single 5 images, after their re-centering on the photometric 
nucleus. With the median average of 5 images, all the possible 
background stars and detector defects (hot or dead pixels) were 
almost completely erased. For the night of July 5-6, this was not 
possible, since only 2 comet and 2 sky images were recorded. In 
this case, the background stars and detector defects were erased 
manually. Although for morphology studies this was acceptable, 
it prevented any precise quantitative measurements for this par- 
ticular night. 

The same procedure was applied to the standards stars. From 
the reduced standard star images, photometric zero points were 
derived for clear nights (using aperture photometry and the pro- 
cedure described in Boehnhardt et al. (2007) for the EMMI im- 
ages). 

3.2. Residual background flux removal 

The presence of a constant residual background was checked and 
corrected by measuring the function SAf at large projected nu- 
cleocentric distances, p. EAf, derived from the Afp introduced 
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Table 1. Filter characteristics and summary log for the comet observation. The table lists on the top the filter name, central 
wavelength, bandpass FWHM and the instrument, as well as, in the bottom, the starting exposure time for the respective filter 
imaging performed before (- hours) or after (+ hours) the impact. The sky conditions are indicated in the last line of the table, 
using the following abbreviations: CLR = clear sky, THN = thin cirrus , COUT = clouded out. On night 5-6/07 the table lists JHK^ 
imaging from the VLT using ISAAC, t-t^) is the time difference between the observation epoch and the DI impact time at the comet. 
It is given in hours and minutes, with negative values before the impact and positive ones after the event. 
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by |A'Heam et al. | (|1984l|, describes the dust albedo (A) multi- 



plied by the total area covered by the solid particles in an annulus 
of radius p and unitary thickness. It is equal to 2npAf, where f 
is the average filling factor of the grains at the projected dis- 
tance p. Note that the definition given here is slightly different 
from the the original one given in Tozzi and Licandro ( 2002| l 
and Tozzi et al. ( 2004| l, even though the physical meaning is the 



same. Assuming a simple outflow pattern, i.e. geometric attenu- 
ation and expansion at constant outflow velocity of the cometary 
dust, the ZAf function should be independent of p. In this case, 
a small residual background (for instance, from incomplete sky 
subtraction) would introduce a linear dependence of the function 
with p. Hence, by applying a trial and error procedure, the resid- 
ual background can be removed from the reduced images such 
that the SAf function becomes constant at large p. This proce- 
dure does not affect the detection of changes in the cometary 
activity, since the latter introduces an "expanding bump" in the 
profile (see section 4.2.2), a very different behavior from the lin- 
ear dependence with p introduced by uncorrected background 
subtraction. This approach for residual background removal is 
still applicable for the observations taken within about a day af- 
ter the impact, since the dust produced by the impact was con- 
fined to cometary distances shorter than 20 000 km and the coma 
flux measured in the SOFI images at larger distances could still 
be used for the above-mentioned calculations. This method of 
the residual background removal is not easily applicable to the 
coma images taken with ISAAC during the night July 5-6, 2005 
since the DI ejecta had already expanded to the edge of the field 
of view. For those observations, we assumed that the integral 
over the position angle (PA) of the function SAf, obtained on 
the opposite side of the ejecta cloud (over PA between and 90 
deg), remained unchanged from night to night. By changing the 
background level in ISAAC images, the flux profiles measured 
in this quadrant before and after DI were forced to be constant 
with p. Indications for the assumption of unchanged appearance 
of the normal coma comes from the analysis of pre-impact ob- 
servations (see [Lara et al.| [2(306) and from the fact that the coma 



signal disappears in the respective PA range when subtracting a 
pre-impact image from a post-impact one (both taken through 
the same filter). 

3.3. Flux calibration 

All the images taken in clear conditions were then calibrated in 
Af using the following formula, derived from A'Hearn et al. 
(1984) 

Af = 5.34 X lO'i (^)' C, X 10-»'*(^"-^') (1) 

where r;, is the heliocentric distance in AU, dx is the detec- 
tor pixel size in arcsec, Cj is the pixel signal in e"/s, Zp and 
Ms are the zero points and the solar magnitude in the used fil- 
ter, respectively. Images taken at non-photometric conditions are 
flux calibrated assuming that the EAf profiles at large p are co- 
incident with those of the day before and/or day after. This was 
justified by the fact that, due to the low dust expansion velocity, 
any change in the dust production would not affect regions at p 
larger than 20 000-30 000 km in 24 hours. Moreover, the coma 
analysis by Boehnhardt et al. |p007| l suggests that no significant 
changes in the flux distribution (except for the DI ejecta cloud) 
took place between July 3 and 10, 2005. The relative calibra- 
tion from consecutive good nights was checked by comparing 
the EAf values of the comet at large nucleocentric distances. 

By a careful examination of all comet exposures, we noticed 
that the images recorded before and 90 hours after the impact 
were very similar and no or negligible traces of DI ejecta were 
found. Hence, in order to increase the signal-to-noise (SN) ratio, 
we computed images (one per filter) of the "undisturbed" comet 
(hereafter called the "quiet" comet) as the median average of the 
comet images taken before and after 90 hours from the impact. 
The standard deviation of the median average is within 2 - 6 % 
for the most part of the comet, i.e. , for regions at nucleocentric 
distances between 2 000 and 50 000 km. In regions closer to the 
nucleus, this standard deviation increases slightly because of the 
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Table 2. Measured Afp for the "quiet" comet (see text) Table 3. XAfo and SAfi best fit results (see text) 



Filter 


Afp(cm) 


Band 


EAfo(cm) 


5:Af,(cm) 




111 ± 11 




345.3 ± 1.3 


10.7 ±4.1 




125 ± 12 




377.9 ± 1.3 


62.8 ± 2.8 




191 ± 19 


R, 


603.5 ± 0.6 


2.5 ±2.8 


J 


228 ± 23 


J 


682.0 ±0.6 


251.0 ± 1.9 


H 


253 ± 25 


H 


758.4 ±0.9 


313.8 ±3.1 


K, 


269 ± 27 




792.9 ± 1.6 


392.4 ± 6.3 



effect of different seeing in the various nights. It also increases 
at distances larger than 50 000 km due to the low coma signal. 

The subsequent scientific analysis of the calibrated images 
is based mainly on SAf profiles and Afp, both easily obtained 
by numerical integration of the flux in the comet images in con- 
centric apertures centered on the nucleus. The physical meaning 
of the EAf profiles is described above. Following the original 
definition by j A'Heam et al. | ( |I984| l, Afp is proportional to the 
average comet flux in the aperture multiplied by its equivalent 
cometocentric projected distance p. This function does not de- 
pend on p for constant outflow velocity. Thus, when using filter 
images taken in the dust continuum bands, Afp is a proxy of the 
dust production rate, Qdust, of the cometary nucleus. However, 
due to unknown dust properties such as the dust size distribution 
and the dust albedo A, it is not straightforward to quantify Qdust 
using Afp measurements of comets. 



4. Results 

4.1. The "quiet" comet 

4.1 .1 . Afp and lAf as a measure of cometary dust production 

The SAf and Afp profiles of the "quiet" comet, derived from the 
observations of 9P in the six continuum bands Uc, B^, R^, J, H, 
Ks during the nights on July 2/3, 3/4, 7/8, 8/9 and 9/10, 2005, 
are plotted in Fig.[T] Various pieces of information on the comet 
dust production can be derived from these profiles. The hori- 
zontal profiles at distances beyond 10000 km from the nucleus 
suggest a steady-state level in the dust production that resem- 
bles homogeneous and isotropic dust expansion in the coma at 
a constant speed. From the existence of jet and fan structures in 
the 9P coma and since the radiation pressure modifies the dust 
distribution in the coma, it is clear that these ideal conditions 
are not fulfilled. However, as long as the jets and fans are stable, 
i.e. they don't change the dust production, the Afp and SAf func- 
tions are constant. The solar radiation pressure may introduce a 
linear dependence of these functions with p, but it becomes no- 
ticeable only at large scales. The Afp values as a measure of the 
dust production of the "quiet" comet are determined at projected 
distances larger than 40 000 km. where the radial profiles in Fig. 
[T]have reached constant values. Results are presented in Table[2] 
The error in the Afp measurements is mainly due to the rela- 
tive photometric calibration error, which is estimated to be of the 
order of 10 %. The Afp values given here are slightly higher than 
the value of 1 12 cm given by Schleicher et al. ( 2005]) for obser- 
vations in the green wavelengths (445-526 nm). They are also 
higher than the value of 102 cm, later revised to 99 cm, derived 
from Rosetta/Osiris observations, using the NAC (Near Angle 
Ca mera) broad-band filters (( |KeUer et al~l |2005| l and ( [Keller et| 
|al. | |2007)). However, as already pointed out by Schleicher et al., 
this may be due to the larger phase angle of the spacecraft ob- 
servations (69°) compared to the measurements from the Earth 
(41°). 



4.1 .2. Signatures of dust sublimation: 

It is evident from Fig. [T]that the near-IR SAf profiles are not 
constant. They increase significantly for distances smaller than 
~ 15 000 km, showing also a little spike very close to the pho- 
tometric nucleus. The spike is probably the signature of the nu- 
cleus convolved with the seeing. However, the SN ratio of this 
signature is too low to derive any useful information. Instead, it 
can be evaluated through the spatial resolution imaging of the 
coma using adaptive optics systems, such as those collected dur- 
ing the impact week using the NACO instrument at the VLT 
(not described here). The slow increase of EAf cannot be due to 
dynamical phenomena (for example, increased cometary activ- 
ity), since the near-IR profiles derived from different observing 
nights look very much the same. Instead the 2Af profiles in the 
visible do not show any evident increase at small nucleocentric 
distances. 

Similar EAf profiles have been found in comet C/2000 WMi 
(LINEAR) (hereafter WMj) ( |Tozzi et al. | |2004l l. At fliat time 
this phenomenon was interpreted as a result of sublimation of 
two kinds of organic grains: one with a lifetime of a!l.3 h and 
the other of » 17 h. Following the analysis of the WMi data, 
the near-IR SAf profiles of 9P for p > 1 000 km were fit by a 
function of the kind 



SAf(p) = SAfo + lAf 1 X e 



(2) 



which contains a constant term EAfo representing the non- 
sublimating (permanent) dust component, and just one decay- 
ing term EAfi representing sublimating grains and character- 
ized by the length-scale Li. The fit achieved for 9P is very 
good and shows length-scales Li of similar value for all three 
near-IR bands, i.e. 6300 + 160 km. Using the fixed length-scale 
6 300 km, the fitting procedure was done one more time, now for 
all profiles including those derived from visible data. The best 
fit parameters EAfo and SAfi, including standard deviations, are 
listed in Table [3] Again, the fit gives very good results, with the 
exceptions of 11^ and R^. filters, where SAf i has values close to 
zero. 

Interesting trends appear when plotting the wavelength de- 
pendence of the decaying (SAfj) and the constant (SAfo) terms 
of the fits (see Fig. |2]i. In the near-IR the constant term varies 
only by a factor of 1.16 going from J to Ks, while the decaying 
term changes by a factor of 1.55. This finding may indicate that 
the two solid components are of a very different nature: one (the 
permanent one) composed of refractory grains and the other (the 
decaying one) made of sublimating grains (or dust covered by 
sublimating material) that scatter very efficiently in the near-IR, 
but are inefficient scatterers in the visible light. As for the comet 
WM], taking into account that the length-scale for the density 
was about 25 % longer than that for the column density, and 
assuming an outflow velocity of about 0.2 km s"', the lifetime 
of the sublimating material was of the order of 40000 s (about 
1 1 hours). Assuming that the sublimation is driven directly by 
the radiation, the scalelength and lifetime scale as the square 
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Fig. 1. Afp and SAf profiles determined from the "quiet" comet (see text): f/c (solid line), (long dashed line), (short dashed 
line), J (dotted line), H (Long dashed - dotted line) and (short dashed - dotted line) 




300 



800 1300 1800 
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Fig. 2. Wavelength dependence of the permanent (solid line) and 
sublimating (dashed line) dust components in the coma of comet 
9P before the impact. Error bars of the fit parameters SAfo and 
SAfi in eq.|2]are also plotted, but are too small to be easily seen 
at the plotted data points. 



of the heUocentric distance, r/,. Then the density length-scale 
and the lifetime at r,, = 1 AU should be 3 500 km and 17 800 s 
5 h). That lifetime differs from those found for the volatile 
dust grains in WMi, which, assuming as well an outflow velocity 
equal to 0.2 km s"', were 61000s (* 17 h) and 4700 s (a; 1.3h) 
at the heliocentric distance of that comet (1.2 AU). They scaled 
to 42000 s (a; 12 h) and 3300 s (54 min) at 1 AU. The gi-ain subli- 
mation may not depend directly on the solar irradiation, but may 
depend indirectly on it, through the grain temperature. In this 
case the lifetime and scalelength scale in a more complicated 
way with r/,, dependending on grain size and thermal degrada- 



tion of the element (see, e.g., Cottin et al. 2004 1. 



The spectral scattering properties of the sublimating grains 
in WMi are different from those found here: both WMi com- 
ponents have a good scattering efficiency in R and /, and the 
long lifetime component has a scattering efficiency similar to 
the normal dust. The sublimating component in 9P almost does 
not scatter in the visible. This means that, even if the phenomena 
may be similar for the two comets, the nature of the grains must 
be different. 

In order to investigate the origin of the sublimating dust com- 
ponent, we have tried to determine its 2D distribution in the 
coma by subtracting a radially symmetric artificial coma image 
from each near-IR filter image. The artificial coma image was 
computed using the parameters for the permanent dust compo- 
nent in the respective fit to the original image. Figure [3] shows 




Fig. 3. Difference image of the "quiet" comet minus a synthetic 
one, as computed from a filter exposure of the comet on July 
2-3, 2005 following the procedure described in the text. North 
is up. East to the left. The Sun direction is indicated. The FoV 
of the image is 40000 x 40000km2 and intensity in Af (0-15 
xlO"^), with a logarithmic look-up-table. The image shows an 
enhanced flux level in the southern coma hemisphere due to the 
sublimating dust and above average activity from isolated re- 
gions on the nucleus (jets&fans). It appears as if the sublimating 
grains are emitted between PA 100 and 200°. 



an example for the difference image in the filter. The non- 
uniform flux distribution in the images suggests that the subli- 
mating dust is more confined in the coma sector defined approx- 
imately by PA = 100° to 200°. Given the indicated PA range, 
the appearance of sublimating dust seems to correlate with the 
nucleocentric surface regions of enhanced activity since various 
jet and fan structures are detectable in the same coma region 
(see 'Lara et al. 2006 and|Boehnhardt et al. |2007|. However, the 
jets and fans are also present in the difference images described 
above since the subtraction of a circular symmetric comet im- 
age does not cancel any asymmetries in the coma flux distribu- 
tion. Thus, it is not possible to fully disentangle the jet and fan 
structures from the sublimating dust component, except that the 
former have much lower intensity levels than those of the fading 
grains. 
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4.2. The ejecta cloud 

4.2.1 . Geometry of the ejecta cloud 

For the study of the various effects produced by the DI event, 
the signal of normal activity was removed from the post-impact 
images of 9P by simple subtraction of the image of the "quiet" 
comet, taken through the same filter This processing should re- 
move the non-impact comet coma without introducing new un- 
wanted features from day-to-day variability, since the normal ac- 
tivity of 9P displayed a rather steady-state appearance. The ex- 
panding cloud of solid particles is clearly noticeable in the visi- 
ble and near-IR images until at least July 6-7, 2005 (three days 
after DI). Fig.|4]shows the ejecta cloud in J band as seen 17:29 
and 20:29 hours after the impact. It can be seen in the figure that 
the cloud is initially expanding into the coma sector between PA 
of 120 to 345°. The time evolution of the cloud expansion can be 



Table 4. Effective cross section of the ejecta grains in km^ vs. 
time t - to after impact for different dust filters. 'Effective cross 
section' is defined in the text. 



characterized by using visible broadband imaging ( Boehnhardt 
[etari|2007l l. 

4.2.2. Ejecta dust production 

Integrating the Af in the images difference over the position an- 
gle range of the initial ejecta cloud (PA=120-345°), we obtain 
the 2Af profile of the cloud vs p. The SAf profiles determined 
in the R^, J, H and K^, filters for the night just after the impact 
are shown in Fig.|5] Note that the different extension of the cloud 
profiles in the figure is due to the different observing epochs dur- 
ing which the cloud was expanding in the field of view. 

By integrating these profiles over p, we can obtain the total 
scattering cross section (SA) of the dust ejecta, SA = j (l,Af)dp, 
i.e., the albedo at the phase angle of the comet multiplied by the 
total geometric grain cross section. SA provides useful informa- 
tion to evaluate the number and the intrinsic color of the parti- 
cles produced by the impact and their evolution with time. For 
the first observations of the first night (July 4-5, 2005) after the 
impact, we measured the SA to be 27.6, 27.3 and 34.6 km^ in 
J, H and respectively. Assuming that the scattering proper- 
ties of the ejecta grains are the same as those of the refractory 
component in the pre-impact coma, it is possible to estimate the 
time interval TI necessary to produce the same amount of dust 
by the normal activity: TI - where v is the mean outflow 

velocity, and SAfo are the values obtained for the normal activ- 
ity. To estimate the order of magnitude, any possible differences 
in the scattering properties of the dust grains were ignored and a 
V = 0.2 km s"' was assumed. Then TI is about 5-6 hours (de- 
pending on the filter). For the lower velocity v = 0.1 kms the 
equivalent duration of normal dust production is doubled. We 
conclude that the amount of dust produced by DI and detectable 
in the near-IR about 17-20 h post-impact time is equal to the 
amount of the dust produced during a few hours (maximum half 
a day) of normal activity of the comet just before the impact. 

Results for the effective scattering cross section SA in km^, 
estimated from the dust filter images available to us, are tabu- 
lated in Table |4] The relative error for different filters depends 
mainly on the respective flux calibration uncertainty that is esti- 
mated to be of the order of 10 %. The relative error for different 
measurements with the same filter does not depend on this cali- 
bration, because during the data reduction the EAf profiles of the 
ejecta cloud images were checked to match the reference pro- 
file of the "quiet" comet at nucleus distances beyond 30 000 km. 
Thus, the relative uncertainty in 2Af for different post-impact 
epochs depends only on the accuracy of the matching of the 
"quiet" comet profiles. As seen in Figs. |4] and |5] a rather ac- 



t - to (hh:mm) 




SA(y) 


SA(H) 


SA(K,) 


17:29 




27.6 


27.3 


34.6 


20:29 




17.9 




33.0 


22:32 


27.3 









curate match is achieved and, hence, the relative error of SAf is 
evaluated to be less that 5 %. The values shown in the table are 
in agreement with the 33 ± 3 km^ obtained in the visible range 
from Rosetta/Osiris observations about 40 minutes after the im- 



pact ( Kiippers et al. 2005 1. 



4.2.3. Ejecta velocities 

Contrary to the natural outbursts in comets, in the case of DI 
event the exact starting time of the formation of the ejecta cloud 
is very well known. Hence, assuming that the major dust pro- 
duction by DI was short-term (as suggested by the fly-by space- 
craft imaging; seel A'Hearn et al. [2005) , the radial profiles of the 
ejecta cloud would also reflect the expansion velocity distribu- 
tion of the cloud particles. The solar radiation pressure applies 
some acceleration ^ to the particles, which is inversely propor- 
tional to the particle radius, a. It has been shown that some time 
between 1 and 2 days after impact the dust grains, as observed 
in the visible, reached the turning point in their motion in the 
sunward direction ( |Boehnhardt et al. ] |2007| l due to the solar ra- 
diation pressure. 

So, depending on the grain size, the radial profiles 2Af(p) of 
the ejecta cloud show a memory of the initial projected veloci- 
ties of the dust after the ejection from the nucleus and possible 
further acceleration in the near-nucleus zone by the ejecta gas 
and/or normal gas release activity. Thus, 2Af(p) divided by the 
elapse time since the impact results in some distribution function 
ZAf(v) where v is the mean velocity of the dust grain at the re- 
spective projected nucleocentric distance in the cloud. A typical 
mean velocity distribution is shown in Fig.[6| As pointed out by 
Kiippers et al. ( 2005 | l and Jorda et al~| ( [2007j ) the radial profile 



of the eject cloud can be well fitted with a Gaussian function. In 
Fig. |6] the Gaussian function that provides the best fit to the ve- 
locity profile is overplotted. The agreement is very good for all 
the filters. For comparison, a Maxwellian function is also shown. 
Here the agreement with the data is good in the leading part of 
the function, but is poor for small velocities. 

Because of the radiation pressure (see above), the velocity 
distribution is not the same in all directions; the mean velocity 
is lower in the Sun direction than in other directions. To study 
this effect, the cloud images have been divided in three sectors: 
SI, S2, and S3 defined by their position angle ranges 145-204°, 
205-264°and 265-325°, respectively. The velocity profiles were 
then fitted with a Gaussian function. Results are shown in Table 
[s] where the mean velocity V and its FWHM is given for nights 
July 4-5 and 5-6. The table does not indicate the errors of the 
fitting since they are very small, of the order of ±2 kms"' and 
±4km s ' for the first and second night after the impact, respec- 
tively. 

Table [s] indicates the following: (1) V, as measured in the 
near-IR, is independent of the filter, (2) V depends on the sector, 
it is slowest in S3, but shows very similar values for the three 
near-IR filters, (3) V in sector S3 slightly decreases with time, 
while in sector S 1 it increases with time. 
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Fig. 4. Processed J filter images of the ejecta cloud recorded about 17 and 20 hours after the impact. The "quiet" comet is removed 
from the post-impact exposures. The intensity scale is in Af (0-3 xlO"^), with a logarithmic look-up-table. North is up, East to the 
left, the field of view (FOV) at the comet is 80 000 X 80 000 km^, the comet nucleus position is in the center of the sub-panels. The 
Sun direction is indicated by the arrow. The cloud 20 h after the impact is less bright than 3 hours earlier. The "hole" at the center is 
an artifact from the subtraction processing of the two exposures pre- and post-impact which had different seeing. Seeing variations 
do not affect the cloud intensity outside of the typical seeing radius of 1-2 arcsec, i.e. about 700 - 1 400 km of the nucleocentric 
distance. 



E 
o 

5 




10000 15000 
p(Km) 



20000 25000 




10000 15000 
p(Km) 



20000 25000 



Fig. 5. SAf profiles of the ejecta cloud vs p measured on night 4-5 July 2005 as a function of p. Left panel: / (solid line), H (dotted 
line) and (long dash line) profiles from observations taken 17:29 h after the impact. Right panel: J (solid line), (long dash 
line) and Rq (dotted fine) ^ 20:29 h (22:32 h for R^) post-impact. Note the changes in the J filter profiles between 17:29 and 20:29 
h post-impact time. 



Table 5. Average projected velocities V and FWHM, in m.s- ' , as Vs2 = 96±7, Vs-i = 93±10ms-i for the night July 4-5 and Vs\ 
determined in the three sectors SI, S2 and S3 of the dust ejecta = 130+7, Vsi - 105+9, Vs3 = 89 + 1 1 m s ' for the subsequent 
cloud (see text) for different filters. night. Using simple physical considerations, the projected dis- 

tance covered by a grain in the S3 sector (^Sun direction) is 

t^o Sector SI \ Sector S2 \ Sector S3 553(f) = V,jt - xt^ and, in the sector SI, (^60° from the 

Band (hh:mm) V FWHM V FWHM V FWHM , ^ t/ , . i dv r(:n\^,2 x\r^\, ■ 1 1 

^ ' a ntisun direction^) S 5 1 (t) - yey/+ 5 ^cos(60)xr. With smiple al- 

gebraic operations the Nej results are {Vs-i{t) + 2Vs 1 (0)/3. Using 
the average values given above, the results are Vej = 1 13 ±16 and 
116±16ms', for the nights July 4-5 and 5-6, respectively. The 
typical FWHM is of the order of 75 m s"' for both nights. In the 
visible the eject velocity is estimated to be of the order of 120 
m s Since only a single observation of this kind is available to 
us, we adopt an uncertainty similar to those in the near-lR. 



4.2.4. Effective scattering cross section 

This picture can be explained by the solar radiation pressure, 
since the projected Sun direction at PA ^ 290° falls almost ex- A surprising finding is the rapid decrease in the effective / scat- 
actly in the middle of sector S3. The average values of the mean tering cross section SA of the ejecta cloud by 35 % during the 
velocities, as in the near-IR, are the following: Vsi = 123+12, time interval from 17:29 to 20:29 hours after the impact as seen 
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Fig. 6. The SAf(v) profile, which is proportional to the mean 
projected velocity distribution of the ejecta dust, as obtained 
from J filter imaging during night 4-5 July 2005. The measured 
profile (solid line) is compared to Gaussian (dashed line) and 
Maxwellian distributions (dotted lines). 



in Table|4] This is also noticeable in Fig.|4] which demonstrates 
that the cloud at 20:29-hours post-impact is significantly fainter 
than three hours earlier (17:29 h post-impact) even though it has 
almost the same extension. Also, the change in the radial profiles 
SAf of the J filter data 17:29 and 20:29 h after the impact shown 
in Fig. |5] confirms the change in the scattering of the ejecta par- 
ticles. This cannot be due to calibration problems, because SA 
is computed using the diff'erence between the impact images mi- 
nus the "quiet" comet: calibration problems would give values 
different from zero in regions far from the cloud. However, this 
is not the case, as can be seen in Figs. |4] and |5] 

The cloud color 17:29h after the impact is almost 'gray' from 
J to H, but its scattering efficiency increases by 25 % in K^. 
Three hours later (20:29 h post-impact), the scattering efficiency 
increases by 84 % between J and (observations in H were not 
taken because of the lack of time). Thus, the reduction in the ef- 
fective scattering cross section of the dust cloud in the near-IR 
between 17:29 and 20:29 h after the DI event was accompanied 
by a strong reddening of the dust. 

To study this, the scattering efficiency for the three near-IR 
filters as a function of particle radius, a, was modeled. In the 
model, spherical particles were assumed with a power-law size 
distribution of power -3.1 as derived by Jorda et al. ( 2007) us- 
ing data obtained 40 min after the impact. The refractive index 
was set to 1.65+10.062, that represents a mixture of silicates and 
organics typical for comets ( Jessberge r et al. | ( 1988)1. Results are 
shown in Fig.|7] the left panel gives the scattering efficiency for / 
and Ks bands as a function of the particle radius, while the right 
panel gives the cumulative scattering efficiency of the particles, 
normalized to 1 at large sizes (5 //m). From the figures it is pos- 
sible to see that the particles with radius less that 0. 1 ;um give a 
negligible contribution to the scattering and that 80 % of the total 
scattering is reached for the particles with a < 1 and < 1 .5 yum for 
J and K^, respectively. Note that the particles with a < 0.92 fj.m 
are mainly responsible for the greatest difference between the 
cumulative scattering in J and K^. Indeed, these particles pro- 
vide 77 % of the total scattering in J and only 44 % in /Tj. From 
these results it is possible to derive an important conclusion: if 
the destruction or sublimation of particles with a < 0.92 fim is re- 
sponsible for the above-mentioned decrease of SA in the J band, 
a destruction of 50 % of such particles should decrease the SA 
by 35 %. However, since these particles (a < 0.92//m) contribute 
also to the scattering in (44 %), this destruction should also 



produce a 22 % decrease of the SA in /Tsfilter. Instead, the mea- 
sured decrease for this filter is only 5 %. 

A possible conclusion is that the observed particles do not 
obey a power law size distribution. The results obtained by the 
DI spacecraft spectrometer indicates that the original DI size dis- 
tribution was dominated by particles of a few microns in size 
( A'Hearn et al. |2005| and supporting on-Une material of Lisse et 
|al. ||2006| ). The difference between Uo rda et al. ipOOT] ) and this 
result may be due to a poor sensitivity of the Jorda et al. mea- 
surements to particles of size larger than 1 fim as is expected for 
the measurements in the visual, as pointed out in Jorda et al.. 
Thus, the particles of size a x I fim could dominate in the ejecta 
cloud 17:29 h after the impact and produce a significant contri- 
bution to the radiation measured in J filter It is well known that 
small grains sublimate faster than large ones because they warm 
up more efficiently than the larger ones, since the absorption of 
the solar radiation is proportional to their area and the heating 
is inverse proportional to their volume. This results in the 1/a 
law (see, e.g. |Lamy~ 1974| l), i.e. the smaller the particles are the 
faster they sublimate. This means that 1 fim particles may subli- 
mate faster and be eliminated from the size distribution more ef- 
ficiently than larger particles. Thus, 20:29 h after the impact the 
maximum of the size distribution shifts to larger particles. This 
leaves the / band without the most efficient contributors whereas 
the situation the band remains almost unchanged. This sce- 
nario, consistent with the in situ data, gives a hope that careful 
simulations of the sublimation of Deep Impact ejecta particles 
may provide some information about the sublimation rate of the 
ejecta volatiles, and, thus, may help to identify it. 

An alternative scenario would be a change in the particle 
composition which significantly modified the value of the re- 
fractive index. However, this hypothesis cannot explain the ob- 
servations. First, the most dramatic changes would be expected 
shortly after the impact when the most volatile components of 
the dust, e.g. ice, sublimate. For the less volatile components the 
change should be slow and it is hard to imagine how such a dra- 
matic change happened between 17:29 and 20:29 hours after the 
impact. Second, we are not aware of any material that may be ex- 
pected in comets that has a so significant difference in its optical 
properties between J and filters. 

4.2.5. Color gradient 

Due to the different orientation of the grain velocity vector to the 
radiation pressure force, radiation pressure has different effects 
on the dust grains in sectors SI, S2, and S3 of the ejecta cloud. 
This results in different sorting of dust particles by size for dif- 
ferent sectors which may appear as a difference in the dust col- 
ors. From the SAf(p) profiles the color distribution of the cloud 
is computed and the results are averaged over the PA range of 
the three considered sectors. Since only the near-IR images had 
sufficient signal-to-noise ratio within the nucleocentric distance 
range of 2 000-9 000 km, and only J and observations were 
performed during the first night after impact, only the J - 
color (C) is computed. The results, in percent per nm and per 
1 000 km, for the night July 4-5, 2005 show that the ejecta cloud 
becomes 'bluer' with larger distance from the nucleus. From the 
images - 17:29 h after the impact we have found the following: 
Csi = (4.1 ± 0.1) - (0.31 ± 0.02)%nOT-i(1000km)-i 
Cs2 = (3.8 ± 0.1) - (0.36 ± 0.02)%«m-i(1000km)-i 
C53 = (5.0 ± 0.1) - (0.51 ± 0.02)%«m-i(1000km)-i 
Thus, the grains closer to the nucleus scattered about 40 % 
(4 % X the difference of the central wavelength of Ks and J) 
more efficiently in Ks than in J. However, they scatter with al- 
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Fig. 7. Differential (left panel) and relative cumulative (right panel) scattering efficiency of a cloud of particles, as a function of 
their radius (see text), as observed in the filter J (solid line) and in (dotted line). 



most the same efficiency in sector S3, at about 10000km from 
the nucleus. The strong difference in the gradients between the 
three sectors can be explained by the effect of the solar radiation 
pressure mentioned above. 

Three hours later (i.e. 20:29 h after the impact) we find: 
Csi = ( 7.5 + 0.3) - (0.28 + 0.06)%nm-i(1000km)-i 
Cs2 = (10.7 ± 0.3) - (1.06 ± 0.02)%«m-i(1000km)-i 
Cs3 = (15.5 ± 0.3) - (1.76 ± 0.02)%«m-i(1000km)-i 
As is expected from decrease in the effective scattering cross 
section in the J filter, a significant increase in the dust reddening 
took place close to the nucleus. The scattering efficiency of the 
grains close to the nucleus became a factor 2-3 higher in than 
in J. We also notice a strong change in the spatial gradients of 
the dust reddening: while the gradient in SI changes only by 
about 25 %, in the two other sectors the change is significantly 
more evident, e.g. by factor of 3 for S3. Moreover, the increase 
is non-linear in the inner part of the cloud, as can be seen in Fig. 

m 

These facts cannot be explained by the solar radiation pres- 
sure only. Also, from our estimation of the speed of particles, the 
three hours between our observations cannot bring a significant 
number of new particles: the particles can travel only k 20 km 
farther from the nucleus. Thus, 17:29 hours after the impact and 
3 hours later we observed particles of the same origin. However, 
their properties and/or size distribution had changed dramati- 
cally. This is confirmed by the above change in the scattering 
cross sections. The changes in color likely also resulted from the 
sublimation of particles which had some specific, non-power- 
law size distribution. This dominance of particles of different 
sizes is what determines the color and its gradient. The differ- 
ence in the color gradients may indicate the efficiency of sub- 
limation whose rate increases as particle size decreases (see 
Section 4.2.4). 

5. Discussion 

An important result of these observations has been the discov- 
ery of the sublimating component in the coma of the so-called 
"quiet" comet. This kind of fading grain has been found in comet 
C/2000 WM], where two components, one with a short and the 
other with a long lifetime, were found when the comet was at 
1 .2 AU from the Sun. These components were interpreted as or- 
ganic grains or refractory grains embedded in organic matter 
that sublimated while heated by the solar radiation. The subli- 
mating component discovered in comet 9P seems to be differ- 
ent from both sublimating components found in comet C/2000 
WMi. It scatters very efficiently in the near-IR, but it does not 
scatter at all in the visible. Scaled to 1 AU, its length-scale is 
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Fig. 8. - J color reddening in % nm ' of the ejecta cloud 
in sector S3 as a function of the projected distance from the nu- 
cleus. Lower graph is obtained from observations taken 17:29 
hours after the impact, the upper graph applies to 20:29 hours 
after DI. Dotted and dashed lines show linear fits to the measure- 
ments. While the color remained almost the same at distances 
p = 8 000 - 10 000 km in the cloud, at smaller distances both the 
color and its spatial gradient increased significantly during the 
3 hours between the two observation sequences. Note also the 
non-linear increase of the reddening function for projected dis- 
tances below 3 000 km, as seen in the data obtained 20:29 hours 
after impact. 



of the order of 3 500 km and, assuming an outflow velo city of 
.2km s~', the lifetime is 18 000s (5h). Note that .Cottin et al. 
j] ( 2004 | l suggested that refractory organic grains, namely poly- 
oxymethylene (POM), may be responsible for the distributed 
source of formaldehyde, observed in several comets. They com- 
puted the scalength and lifetime of the POM grains assuming 
photolysis by solar radiation and thermal sublimation. With a 
temperature of grains of 350 K, they computed a POM length- 
scale of the same order of magnitude as that measured in 9P, de- 
pending very little on grain size. It varied from 3 300 to 7 100 km 
for sizes going from 0.1 to lOjum. 

The following results are obtained for the ejecta cloud pro- 
duced by the impact : 



The total amount of the dust, multiplied by the albedo, covers 
a surface of about 30 km^ about 17 h after the impact, but it 
drops dramatically for the J band 3 h later; 
The velocity distribution of the solid components had a 
Gaussian distribution with an average ejected velocity equal 
to 115±16 ms^'with aFWHM of the order of 75ms"'; 
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- The velocities in the projected direction of the Sun are 
smaller than those in other directions and have similar values 
for the near-IR bands. Those in the visible have larger values 
than the near-IR ones; 

- From the observations 17:29 h hours after the impact, the 
near-IR color of the grains close to the nucleus was found 
to be very red and shows a strong gradient with the nucleo- 
centric distance with the highest values in the Sun direction. 
Three hours later, the near-IR color becomes even more red 
and the gradient with p in the Sun direction increases by a 
factor ^ 3. 



The Gaussian velocity distribution of the particles is puz- 
zling. In the case of gas drag produced by an explosive event 
that has a timescale shorter than the acceleration time scale, the 
velocity of a particle with radius a should change as a '. This, 
combined with a power-law size distribution of power equal to 
-3.1 ( Jorda et al. | 2007 ) would produce a velocity distribution 
very different from a Gaussian one. If the lifetime of the explo- 
sive event was longer than the timescale of the acceleration the 
ejection velocity distribution with particle sizes would follow a 
more complex law, but always size dependent (see, e.g. |Gombosi 



et al. ||1986) ). In the case of a natural outburst the velocity dis- 



tribution has been found to follow a Maxwellian one (see, e.g., 
[Schulz et al. [ |2000| ). On the other hand, a thermal acceleration 
of the dust, that would give a Maxwellian distribution, is ex- 
cluded because it would require a temperature excessively high 
to accelerate a grain with the mass of the order of 10"'^ g at a 
velocity of a hundred m s"' . Models of ejecta with grains with a 
power law size distribution, accelerated by gas drag, give a ve- 
locity distribution far from a Gaussian one. Only in the case of 
grains with almost the same size would the velocity distribution 
follow that of the grains and it can also became a Gaussian one. 
The average projected velocity found her e is in good agreement 



Feldman et al. 



( |2007| l 



with the value of 11 5 m s" obtained by by 1 
from HST observations in the visible jus t after the impact. It is 
smaller than those noted byjMeech et al. |(|2005) (200+20 m s"') 



and Schleich er et al. ( 2005 i (220 ms"') who refer to the veloc- 
ity of the leading part of the cloud. For the latter we find 250-300 
ms"', in fairly good agreement with their values. 

However, our result for the Gaussian velocity distribution of 
the expanding ejecta cloud differs from that of Jorda et al. ( 2007| l 
(190ms ' with FWHM - 150ms '), derived from Osiris obser- 
vations on-board Rosetta. This difference cannot be explained 
with the different viewing geometry of Rosetta spacecraft with 
respect to an observer on Earth (the difference in aspect angle 
is just 20°). The Osiris and our ESO measurements, obtained at 
similar wavelengths, are in a good agreement for the total light 
scattering area of the dust ejecta (33 km^ by Osiris and 27.3 km^ 
for our observations). An explanation for the different veloc- 
ity of small and faster grains resulting from dust sublimation a 
few hours after the impact is thus very unlikely. 

The dramatic changes in the total amount of dust observed 
in J band is another very puzzling result. It is associated with a 
strong change in the color of the cloud and its gradient with p. A 
possible scenario would be a sublimation of the dust grains con- 
taining slow-sublimating volatiles such as organics. Water ice is 
excluded because its sublimation timescale is much shorter than 
17 h (Hanner 1981[ l. If the particles of a few micron size dom- 
inated in the original size distribution (as was found from the 
in-situ data) then sublimation of such particles would result in 
elimination of the most abundant 1 micron size particles, which 
are the most efficient in the J band. This will manifest in a signif- 



icant decrease of the brightness in J band as well as by a change 
in the dust color. 

This means that the scattering area SA of the grains just af- 
ter the impact must have been much larger than that measured 
at 17:29 h post-impact, i.e., the quantity of the solid component 
released by the impact may have been of an order of magni- 
tude larger than the 5-lOh of normal activity derived in Section 
4.2 without considering sublimation effects. However this is not 
confirmed by the results of Kiippers et al. ( 2005 i who found a 
light scattering cross section of the cloud in agreement with that 
found here (see above). It is not possible to obtain more infor- 
mation from the results presented here, but about 17 hours after 
the impact the size distribution of the grains seems not to follow 
a power law (if it ever did), but seems "monochromatic". 

There is no obvious correlation between the organic grains 
found in the 'quiet' comet and those supposed to be present in 
the ejecta cloud. 

6. Conclusions 

From observations of gas emission-free regions of the comet 
9P/Tempel 1 made before and after the Deep Impact event, the 
scattering characteristics and the velocity of the ejecta cloud pro- 
duced by the impact have been measured. 

Seventeen and a half hours after the impact, the total area 
covered by the grains of the ejecta cloud, multiplied by their 
albedo, was 27 - 35 km^ in JHKs . Three hours later, it dropped 
to ^ 18 km^ in J, but remained almost constant in K^. During this 
interval of time, the J - color gradient with the nucleocentric 
distance p also changed significantly, increasing by a factor of 
three in the direction of the Sun. 

The projected average velocity of the ejected cloud measured 
in the near-IR was 1 15 ± 16 m s ' , and was found to be indepen- 
dent of the filter used for the observations and the position angle. 
Its distribution was very similar to a Gaussian. 

It has been shown that all these results cannot be explained 
assuming that grain size has a power distribution nor it is possi- 
ble to assume that grains are ejected by gas drag. 

While the mechanism of grain ejection is difficult to explain, 
the behavior afterward can be justified only with the presence 
of a; 1 yum size organic grains that are sublimated by the solar 
radiation. 

From the pre-impact and late post-impact observations, the 
presence of a sublimating component has been detected and in- 
terpreted in terms of organic grains that sublimate because of the 
solar radiation. 
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